Pair Density Wave in the Pseudogap State of High Temperature Superconductors 
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Recent scanning tunneling microscopy (STM) experiments of Bi2Sr2CaCu20g+5 have shown ev- 
idence of real-space organization of electronic states at low energies in the pseudogap stateQ. We 
argue based on symmetry considerations as well as model calculations that the experimentally ob- 
served modulations are due to a density wave of d-wave Cooper-pairs without global phase coherence. 
We show that STM measurements can distinguish a pair-density-wave from more typical electronic 
modulations such as those due to charge density wave ordering or scattering from an onsite periodic 
potential. 
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One of the major challenges in condensed matter 
physics today is to understand the electronic phase di- 
agram of high temperature cuprate superconductors. At 
the center of the debate is the nature of electronic states 
in lightly hole-doped cuprates, which has attracted a 
great deal of attention because of the observation of 
a pseudogap in the density of states (DOS) in these 
compounds To understand this phenomenon, some 
have focused on the smooth evolution of the pseudo- 
gap into the superconducting gapQ, the anomalous high 
frequency conductivity^], or the Nernst effect 0, 0, 
It was proposed that superconducting pairingfluctua- 
tions persist above Tc in underdoped samplesQ, 0, 
and are responsible for the strange electronic behavior in 
this regime. Others have argued that in the underdoped 
regime the electronic system has a strong propensity for 
ordering hence proposing the possibility for antiferro- 
magnetic order [llj, combined antiferromagnetic and su- 
perconducting fluctuations [l^ . stripe order 
charge density wave id ll^. staggered flux phasejlSj and 
d-density wave(li|. 

Recently, mapping of the electronic states with the 
scanning tunneling microscope (STM) has demonstrated 
a link between the observation of real space electronic 
patterns and the pseudogap in DOS of underdoped 
Bi2Sr2CaCu208+5 samples y. In this experiment, STM 
was used to detect a dispersionless modulation of the 
low energy electronic states, showing an unusual en- 
hancement of the intensity of modulated patterns within 
the pseudogap energy scale. The pseudogap modula- 
tions bears similarity to those observed near the vortex 
cores 20j, but seems to be distinct from the STM features 
in the superconducting stateP, 0,1131, which exhibit en- 
ergy dispersion. The origin of DOS modulation in the su- 
perconducting state has been for the most part attributed 
to quantum interference of quasi-particlesjl^, although 
potential relevance of electronic ordering has also been 
consideredEilElliililEilii. In Refs.Elllii, it was 



argued that the energy dependence of the tunneling in- 
tensity can be used to identify the type of ordering in the 
superconducting state. 

In this Letter, we focus on the STM resolved mod- 
ulations of the DOS in the pseudogap state and pro- 
pose that these are a consequence of a novel density 
wave associated with d-wave Cooper pairs, first proposed 
in the context of the charge ordering in the vicinity of 
vortices 20, 24] . This state, illustrated in Fig. (1) of 
Ref. [2J has a rotationally symmetric charge period- 
icity of 4a X 4a near doping level x = 1/8. This state 
also arises naturally in the plaquette boson approach of 
Altman and Auerbachl29| or in the large N approxima- 
tion of the t — J model'28']. It differs from the stripe 
state in terms of the rotation symmetry, and differs from 
the Wigner crystal state of the holes, which would have 
a charge periodicity of \/8i x \/8a at the same doping 
level. We argue based on both symmetry and model cal- 
culations that a pair-density-wave without global phase 
coherence provides a natural explanation of the energy- 
dependence of experimentally observed patterns in the 
pseudogap state/]. Furthermore, we propose analysis 
of STM data which will allow a pair-density-wave to be 
clearly distinguished from the typical particle-hole charge 
density wave (CDW) or other modulation due to electron 
scattering from a periodic potential. This analysis relies 
on the symmetry of the Fourier transform of the tunnel- 
ing density of states, p{Q,E), which we predict to be 
approximately symmetric for a pair-density-wave when 
E —E, while in case of CDW or potential scatter- 
ing p(Q, E) is approximately anti-symmetric. Our nu- 
merical results shows that this signature of pair-density- 
wave (PDW) state survives breaking of the particle-hole 
symmetry and provides a sufficiently clear-cut method to 
test our proposal. Further experimental confirmation of 
PDW scenario would not only establish the presence of 
an unusual ordering phenomenon but it would also vali- 
date the view that the pseudogap state has indeed short 
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range pairing correlations. 

To model the experimental findings, we study quasi- 
particle spectrum based on solutions of the Bogoliubov- 
de Gennes (BdG) equations for a d-wave superconductor, 
while taking into account classical phase fluctuations to 
model the pseudogap statefsof. 
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with HoiJnir) = -tJ2s'>Pn{r + S) - ^(r)?/;„(r) and 
AV'ra(r) = A(r, (5)'(/;„(r + (5) subjected to the self- 
consistent condition 

Air, 5) = -!:l^^[«„(rK(r')+u„(r')<(r)] 



X tanh 



(2) 



where r' — r + S and S denote a nearest-neighbor vector, 
tpn{^) can be Un(r) or Vn{r). /3 — l/fcsT is the inverse of 
temperature T. We assume V{r,r') — {V{r) + V{r'))/2 
and focus on the periodic pairing and chemical potentials, 
namely, V{r) = Vq - J2q cos(Q • r) and /x(r) = 

- Eg ^KQ) cos(Q • r). 

To describe the pseudogap state above Tc, we consider 
only classical phase fluctuations i.e. the phase field 0(r) 
is independent of time. The pairing operator takes the 
form 

A = le"*'^"^/^ [A(r)d(p) + d(p)A(r)] e^'^^''^/^ (3) 

where the differential operators d{p) — 
2[cos(jixa) — cosipytt)]. In the superconducting state, 
the phase field (/>(r) is rigid while in the pseudogap its 
fluctuations destroy long range phase coherence. Note 
that in the superconducting state, even in the presence 
of pairing modulations, the spectrum remains gapless at 
four points on the Fermi surface. 

The local DOS (LDOS) is a functional of the phase 
configuration, (/'(x), given by 



p(r,i?,{0(x)}) = 



\u„{r)\^f'iE^^E) 



+ \vn{r)\^f{E„+E) 



, (4) 



where the (M„(r), w„(r)) and En are respectively eigen- 
functions and eigenenergy of (Q) for the phase configura- 
tion 0(x). f'{E) = -p/Acosh^{l3E) is the derivative of 
the Fermi distribution function. The Fourier transform 
ofp(r,i?,{0(x)}) is 

p(q, E, {0(x)}) ^^Y. • '■)^('"' ^' 



with the number of the sites. 

Before presenting the numerical results, we analyze the 
symmetry of p{q,E, {^(x)}) when the average of chem- 
ical potential /x(r) is zero. First, we consider the ef- 
fect of a CDW described by a periodic chemical poten- 
tial fi{r) oscillating around zero average with periodicity 
L^. The Hamiltonian has the following symmetry: if 
(M„(r), Vn{r)) is an eigenstate with energy En and phase 
configuration (/)(x), then (r + L^/2), m* (r -|- 

L^/2)) is an eigenstate with the same energy £"„ for the 
phase configuration 0'(x) = (/)(x -|- L^/2). We used the 
fact that the factor (^—l^+y introduces a minus sign to 
both the hopping and the pairing terms leaving the chem- 
ical potential term untouched, while the shift by L^/2 
changes the sign of the chemical potential term. In addi- 
tion the standard Bogoliubov-de Gennes symmetry has 
been used. For the superconducting state where 0(x) is 
x-independent one can easily show yo(Q^, E, {6(x)\)^y\\. 



p{Q^,E,{^}) = -piQ^,-E,{^}) 



(6) 



where • = 27r. In a pseudogap state where 
the phase coherence is destroyed, we must carry out 
a weighted average over all of the classical phase con- 
figurations. We assume that to a good approximation 
the probability distribution P({(/)(x)}) is invariant un- 
der translation, i.e. P({(?!)(r)}) = P({0(r L^/2)}) (e.g. 
XY- model). We thus conclude that the symmetry © 
survives in the presence of the classical phase fluctua- 
tions. 

Now in contrast to CDW, we consider the case of pair- 
density- wave. In the case of the periodic pairing potential 
with periodicity and zero uniform chemical potential, 
(— ?;*(r), u* (r)) is an eigenstate with energy En, 
which makes the Fourier transform of the LDOS an even 
function of energy 



p{q,E,{4>{^)})^piq,~E,{^i^)}), 



(7) 



where q is an arbitrary wave vector. 

In the presence of modulation of both chemical po- 
tential and d-wave pairing potential, the above ar- 
guments do not hold in general. However, in the 
special case of the chemical potential oscillating with 
= (27r/8a, 27r/8a) and pairing potential with = 
(27r/4a,0) and (0, 27r/4a), one can show that the phase 
averaged LDOS satisfy p{Qd,E) = p{Qd,—E) while 
p{Qfj,,E) = —p{Qf^,—E). We beheve that this case is 
relevant for the actual experiments in BSCCO. 

The dichotomy of the local and non-local modulation 
serves as our point of departure in analyzing perturba- 
tions which break particle hole symmetry. If the depar- 
ture from the symmetry is not too large, the above ar- 
guments are still useful as the odd or even character of 
/9(Q, E) is true approximately. An immediate check is the 
tunneling density of states at the wavevector correspond- 
ing to the BSCCO superstructure. We expect it to enter 



3 




-40 40 

E{meV) 



FIG. 1: (a) The average DOS. (b) The Fourier components of 
the LDOS at the wavevectors of the pair-density-wave. T^ '^ 
is the mean-field transition temperature. 



as a local periodic scattering potential, making p(Qp, E) 
odd under E — > —E. On the other hand, the charge den- 
sity wave of d-wave bosons (see Fig. 3 of Refjlll) should 
serve as an effective periodic pairing potential that is in- 
herently non-local. Thus p{Q,d,E) should be even under 
E^ -E. 

We shall now turn to the numerical solutions of the 
BdG equations Q for the superconducting state with 
uniform phase. The effect of phase fluctuations will then 
be described later in the text. We set Vq = 0.5/: so which 
for ^(r) = and V{r) — Vq yields a gap A of about 0.32i 
at low temperature. We choose t = Q.12heV such that 
A = AOmeV, which is relevant for slightly underdoped 
BSCCO. In Fig^we plot the results for a pair-density- 
wave with AV{Qd) = 0.2^0 with = (27r/4a, 0) and 
(0,27r/4a). As expected, p{Qd, E) are even functions 
of energy. At low temperature compared to the mean- 
field transition temperature T*^^, p(Qrf, E) shows peaks 
at the energy within the superconducting gap and zero 
crossings at a energies above the gap, which is consis- 
tent with the experimental data of Howald et al.'21*\. 
At higher temperature, T = 0.86r*^^, piQd,E) has 
a peak at zero energy and decreases as the energy in- 
creases. Again, this is consistent with the experiment 
0. In the case of the periodic chemical potential with 
A^(Q^) = 4meV and = (27r/4a,0), p(Q^,£;) plot- 
ted in FigEl^a) as a sohd line. As expected p{Qf^,E) 
is an odd function of energy. As discussed later, this 
feature enables us to exclude the modulation of chemical 
potential as the direct cause of the observed peak around 
(27r/4a,0). 

So far, our calculation and symmetry arguments as- 
sume the particle-hole symmetry. As an illustration of 
the effects of weak symmetry breaking, in Fig|21we dis- 
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FIG. 2: (a) Fourier components of LDOS with periodic chem- 
ical potential and uniform d-wave pairing, (b) Fourier com- 
ponents of LDOS with periodic d-wave pairing and uniform 
chemical potential. The temperature is T = 0.86T^'^ for 
both cases. 
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FIG. 3: The Fourier components of LDOS for a more re- 
alistic model with next-nearest-hopping t' = —0.3t, chemi- 
cal potential average fj, = —l.Ot, the periodic modulation of 
pairing with wavevector (27r/4a, 0) and (0, 27r/4a) and a pe- 
riodic chemical potential with wavevector (27r/10a, 27r/10a). 
The amplitudes of modulations are AV{Qd) = 0.2Vb and 
A/i(Q^) = O.lt. The dash line is the average DOS. The 
solid and dash-dot lines are the Fourier components of the 
wave-vectors of pair-density-wave and CDW, respectively. 



play p(Q, E) for a finite chemical potential p = —20meV. 
In Figl^lwe include finite next-nearest neighbor hopping 
term t' = — 0.3t and a chemical modulation around a 
finite average with a wavevector (27r/10a, 27r/10a) corre- 
sponding a periodicity about 7a in the (tt, tt) direction. 
In both plots, remnants of the expected symmetry are 
clearly visible. In particular, in the case of the periodic 
p{r), p{Qf^,E) has odd number of zero-crossings below 
the gap energy, while in the case of pair-density-wave 
p{Qd, E) has even number of zero-crossings at energies 
close to the gap. 

We now include the effect of the phase fluctuations 
within the model of Franz and Millisj30|, where the 
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FIG. 4: The LDOS in the pseudogap state calculated from 
Eq.® and FigOlfor Qd = (27r/4a, 0) and (0,27r/4a). Here, 
A ~ AOmeV is the pseudogap. The solid dots are the experi- 
mental data of (27r/4.7a,0) peak of Vershinin et aZQ. 



LDOS of a pseudogap state p(r, E) can be calculated 
as an average over the phase fluctuations that are slowly 
varying in space 



p{r,E) = J dr,P{ri)p{v,E-ri), 



(8) 



where /9(r, E) is the LDOS for the superconducting state, 
•q is just the Doppler energy shift due to the presence 
of a uniform supercurrent and P{vi) is the probability 
distribution of 77 that can be calculated within 2D XY 
Hamiltonian to yield [sof 



P{-q) = V27rVF(r)e 



(9) 



with W{T) a temperature-dependent parameter of the 
order of magnitude of the pseudogap. Since P{vi) is trans- 
lationally invariant, the argument regarding the symme- 
try of the LDOS is stiU valid. In FigH we plot the LDOS 
calculated from Eq.JHl) and FiglSlwith different values of 
parameter W . The curve with W = 0.5A agrees with the 
experimental data of Ref. 0] reasonably well. This value 
of W is consistent with the one obtained in Ref. [soj . 

In conclusion, we show that the Fourier transform of 
the tunnelling density of states recently obtained from 
the STM experiments, p{Q,E), contains sufficient infor- 
mation to distinguish the pair-density-wave state from 
the conventional CDW state. Near the doping level of 
X = 1/8, the pair-density- wave has a rotationally in- 
variant periodicity of 4a x 4a, which is close to the ex- 
perimental observation. We are able to reproduce the 
most salient future of the non-dispersive peak as shown 
in Fig.2C of Vershinin et al.^, namely, the intensity of 
the new peak picks up as the pseudogap opens. Our work 
also makes an experimentally testable prediction about 
the behavior of p{Q, E) ioi E < 0, which has the opposite 
behavior for the pair-density-wave state and the conven- 
tional CDW state. The preliminary analysis of the exist- 
ing data has confirmed that the new non-dispersive peak 



is spatially particle-hole symmetric |33, as predicted here. 
The pair-density- wave state therefore offers a unified the- 
oretical explanation of the STM experiments both near 
the vortex coreji^l and in the pseudogap state. This state 
plays an essential part in understanding the zero temper- 
ature global phase diagram of the high Tc cuDrates|33|. 
and our current work shows that it is also the domi- 
nant competing order in the pseudogap state. This state 
could be favored in the pseudogap regime of underdoped 
cuprates because of the combination of the low superfluid 
density and strong pairing, which is consistent with the 
previous theoretical ideasP, |3 ■ 
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